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We derive Bohm's trajectories from Bell's beables for arbitrary bipartite systems composed by dissipative 
noninteracting harmonic oscillators at finite temperature. As an application of our result, we calculate the 
Bohmian trajectories of particles described by a generalized Werner state, comparing the trajectories when the 
sate is either separable or entangled. We show that qualitative differences appear in the trajectories for entangled 
states as compared with those for separable states. 
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Entanglement phenomenon is possibly the most striking feature of Quantum Mechanics, playing a key role in quantum in- 
formation processing and quantum computing. The striking feature of entanglement was experimentally turned possible after 
the seminal paper by Bell [1], and since then efforts from experimental and theoretical research to demonstrate the violation of 
, Bell's inequalities have been undertaken. To unequivocally demonstrate violation of Bell's inequality, a major development in 
£NJ ' experimental techniques has been carried out to produce entangled photons from the cascade atomic process [2] to the paramet- 
ric down-conversion processes [ 3 ] . Massive entangled particles have also been produced through radiation-matter interaction in 
cavity QED [4] and trapped ions [5]. In the latter case, a controlled entanglement of 14 quantum bits has been recently generated 
enabling the implementation of the largest quantum register to date [6]. The violation of a form of Bell's inequality has also been 
verified with massive entangled particles within an ion trap @]. Parallel to the experimental achievements, theoretical physics 
has been struggling with recently advanced striking features of entanglement, such as the derivation of separability criterion for 
density matrices (U, entanglement sudden death [9], and quantum discord fioll . 

Bohmian mechanics is a theory equivalent with orthodox quantum mechanics having the advantage of providing ontological 
, ^ , ', meaning for the quantum particle trajectory lUTll . Bohmian mechanics assumes that the complete description of particle systems 
is provided by its wave function ^ and its configuration Q = (Qi, .., Qn) € M. 3N , where Q a is the position of the a-th particle. 
t-H | While the wave function ^>(Q) evolves according to the Schrdinger's equation, the motion of the particles evolves according to 
the equation m a dQ a /dt = him { 1 i>~ 1 d , i>/dQa}- Since only the average trajectories are experimentally accessible, the particle 
positions are the "hidden variables" of Bohmian mechanics. 

In the present work we focus on entangled states of bipartite dissipative systems from the Bohmian trajectories perspective 
as formulated by Vink's extension of Bell's beables flafoll . Our goal is to verify what happens to the Bohmian trajectories of 
a given bipartite state at the instant occurring separability. To achieve our goal, we first derive Bohm's trajectories from Bell's 
beables for arbitrary bipartite states under thermal reservoirs at finite temperature. On this regard, we note that an extension 
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of Bell's beables that encompasses dissipation and decoherence for one particle state has been advanced Ill4ll . where the dif- 
fusive terms in Nelson's stochastic formalism are naturally incorporated into Bohm's causal dynamics. Summarizing, here we 
generalize the approach of Ref. lfl4ll to include dissipative two-particle states, thus allowing us to study Bohmian trajectories for 
^ . correlated quantum systems described either by pure or mixed states under dissipation at finite temperature. 

C3 . II. BOHM'S TRAJECTORIES FOR TWO-PARTICLES DENSITY MATRICES 

In this section we derive Bohm's trajectories from Bell's beables to encompass two entangled noninteracting particles under 
independent reservoirs at finite temperature, whose master equation is 



dp{t) _i_ 
dt ~ ~h 



[H, p (t)] + ^ -y [( n » + l)( 2 «aP«a - a f a a a p - pa f a a a ) 
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with a a and oj, being the usual annihilation and creation operators in Fock spaces, j a is the corresponding dissipative rate 
with reservoir average thermal photon number n a , and a = {1, 2} refers to each particle of the system and its corresponding 
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reservoir. This master equation generates the probability density P ni n 2 (t) = {fm > X« 2 1 P (t) [fm , X« 2 ) sucn mat 

u dP nin2 _ , , dp jt) = ^ 

Q£ \ < y 3 ni!Xri2l l ( r«l>Xn2/ — / y >^nin 2 mim 2 - vA) 

mi ,m2 

Now, using the Eq. (fl} and completeness relations we obtain 

J ni n 2mi m 2 = 2Im| {f ni ,Xni\H \tp mi ,Xm 2 ) (¥>mi > Xm 2 I /=>(*) |¥>«i,X™ 2 ) } 

+ | X] 7 a (" Q + l)|2{^„ 1 ,Xn 2 l a c«l^m 1 ,Xm 2 )(^m 1 ,Xm 2 |p(t)|^fc 1 ,Xfc2>(^ftl.Xfc2l a ll^™l.Xn 2 > 

- <Vni , Xn 2 | 4 |'y3m 1 ,Xm 2 ) (^m, , Xm 2 I a a |Vfcl)Xfca> {^fcl ! Xk 2 I |¥>«1>X« 2 ) 

- (<Pni , X«2 I P(t) I^DXmj) (Vto! , Xm 2 I oj, |^fc 1 ,Xfc 2 ) (V?fcl , Xfc 2 I IWl,X» 2 ) } 

- (Vni , X«2 I l^mi,Xm 2 ) {<£mi , X™ 2 I a l IVfcl.Xfc 2 ) fa*! ! Xfc 2 I l^«l,Xn 2 ) 

~ (fn! , Xn 2 | /0(t) |^ roi ,Xm 2 ) , Xm 2 | aa |<^fei,Xfc 2 > (V^l > Xfc 2 I a i IVni)Xn 2 > }• (3) 

As the classical counterpart to the continuity Eq. d2j we write the following master equation for two particles 

d p »i"2 - V fT P -T P ) (4) 

/ J \- L n\n 2 7n-im 2 1 m\rn 2 ± m\m 2 n\n 2 ± n\n 2 ) j V^v 
mi ,m 2 

where T nin2mim2 dt is the transition probability governing jumps from states \<p ni ) and |xmi) to \(p n2 ) and |xm 2 )> respectively. 
The quantum and stochastic formalism meet a common ground through the mixed quantum-classical equation 

Jn 1 n 2 mim 2 _ rp p rp p ,r\ 

. 1- n\n 2 m\m 2 im\m 2 -L m\m 2 n\n 2 in\n 2 \P) 



which admits the particular simplified solution 



Jnin 2 mim 2 ^ 



T — i nP m , m ~ ' "«i" 2 mim 2 ff -\ 

1 U , Jn 1 n 2 m\m 2 _ U 

Next, we assume that the entangled systems are non-interacting harmonic oscillators of frequencies oj a and masses M a , each 
modeled by the Hamiltonian 



H M _ Pa , M «" a r 2 (7) 



p Q being the canonically conjugate momentum to the coordinate variable x a , and the total Hamiltonian is H = if' 1 ) + ff( 2 '. 
Additionally, we consider both entangled systems to be described by the arbitrary general mixed state 

pW-E^I^WXV^WI, (8) 

u 

with (¥>„i, X„ 2 \ip u (t)} = d,„ 2 - 

In Vink's extension of Bell's beables ITOl . where all the degrees of freedom must be discrete and finite, the position is 
restricted to sites of a lattice which, in the one-dimensional case, becomes being integers and e is the lattice 

distance. To extend Vink's approach to two particles (the continuous limit is recovered taking e — > 0) we must z) write the 
smooth wave functions in the coordinate representations as (\<p na ) = l^na)) V'ni n 2 = -^ni n 2 ex P [l'-'ni »J > wnere V'ni n 2 = 
V>" (ar m , x„ 2 , i), i?^ „ 2 ee i?" (x ni , x„ 2 , t), and n2 ee S u (x ni , x„ 2 , t); expand t/; 11 to first order in e, i.e., 

Ci±l ) n 2 =Ci,n 2 ±£A 1 d,n 2 (9) 

<,„ 2 ±i=Ci,n 2 ±eA 2 C iin2 (10) 
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with 
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z'z'j) substitute a na „ la — [M a aj a x ma S na ,m a + fr(8 no ,+i,m a — &n a ,m a ) / z]/ (2hM a uj a ) 1 ^ 2 m the transition matrix given by 
Jn 1 n 2 m 1 m 2 /K neglect terms of order 0(e 2 ) and higher (taking e a = e for simplicity); z'v) take the limit e — > with 
x m2 ±i — em 2 ± e. After a straightforward but length calculations we obtain (a, j3 = 1, 2): 



Jnir 



hj a (2n a + 1) 1 
2 A4 a to a R u 



^m a —l,n a ) Urn/j ,n 



)S„ 



Next, by defining 



c a (t + dt) ~ x a (t) + e[(k — m)d at i + (I — n)8 a ^ 



(12) 



(13) 



with (k—m) — J2mn Tmkni (k — m) and (I — n) = ^2, mn T m knl {l—ri)dt, such that for the forward movement n\ — m\ = 1 
i.e., k > mil > n) for particle 1 (2), Eq. ( fT3l results in 

dx a (t) 1 



dt 



J2 u Pu(R u (xi,x 2 ,t)y 



(x 1 ,x 2 ,t)) 



1 dS u {x u x 2 ,t) h la {2n a + \) 



M a dx a 



dR u {x 1 ,x 2 ,t) 



2M a uj a \R u (xi,x 2 ,t) dx a 



The above equation can be rewritten more compactly in the following way, using p [x\, x 2 , x[, x' 2 ,t) — (xi, x 2 \ p (t) \x[, x' 2 ), 



dx a (t) h 
~dt ~ Ma 



Im [d Xa p (xi ,x 2 ,x' 1 ,x 2l t)} 
p(x 1 ,x 2 ,x' 1 ,x' 2 ,t) 



x\—x l 
x 2 =x' 2 



h^ a (2n a + 1) 
2M a uj a 



Re [d Xa p (xi , x 2 , x[ , x' 2 , t)] 



p(xi,x 2 ,x[,x' 2 ,t) 



X2=x' 2 



(14) 



We note that our main result, Eq. (TBI , generalizes that one obtained by a different approach in Ref . 111511 for density matrix 
of individual systems without dissipation (7^ = 0). In the following, we shall use Eq. (TBI) to calculate trajectories of quantum 
particles when the entanglement either is present or absent in the joint state. 



III. GENERALIZED WERNER STATES AND MOTION EQUATIONS 

With the Bohmian equation of motion for two particles in hands, Eq. (fl4] i. we next assume that the entangled state is prepared 
in the generalized Werner state 

p = e|^ ± ) 12 <^ ± | (15) 

where I = Ii <Ei l 2 stands for the identity operator and 

^ ± > 12 =a|00)±6|ll>, (16) 

with a and b being complex constants (|a| 2 + \b\ 2 = 1). 

To compute Bohmian trajectories we must find the solution p(t) in the presence of losses due to a thermal reservoir. We then 
will specialize to the case of losses at zero temperature, using the method of phenomenological operator approach, as developed 
in Ref. [19], where we define Qy a = 7) 

|0)|0) fl |0)|0) fl 

|l)|0) fl -> e-** |1) |0)^. + x^^^|0) 

If we now take into account that the Werner-like state is uncoupled from the reservoir at t = 0, then using the rules given above 
we can write, after tracing out the reservoir variables, 
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eb 2 



-ft 



1-e 
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1 -e 



5& 2 



1-e 



sb 2 



e" 2 ^* |11) <11| 
1 - el 



(l- e -7*) 2 + i_I(i_ e -7*)^|oo) (00| 



'1 - e-^) 



1 - e 



(|10)(10| + |01) (01|) 



±ee" 7t [ab* e 2lut |00) (11| + ba*e~ 2lut |11) 



(17) 



where we have chosen u; Q = uj. We now assume that both entangled particles, originally represented in abstract Fock spaces, 
are harmonic oscillators within the subspace {0, 1} of the ground and first excited states. This assumption enables us to an- 
alyze, through both particles' trajectories, how entanglement dynamics affects Bohmian trajectories. Considering the scaled 
dimensionless variables x a = yjui/hx a , we write the state (fTTI i in the coordinate representation to obtain p (x%, x 2 , x[, x' 2 , t) as 



p(x 1 ,x 2 ,x' l ,x' 2 ,t) = {^A | 



1-e 



dr 



2e~ 7 * 



ea - 
1-e 



eb 2 



1 - e 



4 
eb 2 



e 2lt x\X\x 2 x 2 



1 - e 



fl - e" 7t ) + ^-r-^ \ \x 1 xi'+~x 2 x 2 n 



Using Eq. ( TT4| ) we obtain, after a straightforward calculation 



where 



dx 2 
dt 



7 



£ (fi,f 2 ;t) 

7 

5 (£i,£ 2 ;£) 



=F— ^4(t) sin(2wt)£ 2 + B{t)x x x 2 + C(t)fi ± Ait) cos{2uit)x 2 
7 

=p— A(t) sin(2u;t)xi + B{t)xi 2 x 2 + C(t)x 2 ± A(t) cos(2ujt)x 1 
7 



2 ' 
2 ' 



(18) 

(19) 
(20) 



G (xi,x 2 ;t) 



1 - e 



ea 4 

I - e 
2 

2e~ 7 * 



e& 2 



1 - e 



(l-e- 7t )' 



l-e" 7 *) +4 
1-e' 



cb 2 



1 -e 



e" 274 ^ 2 ^ 2 



d/ 



(1 - e" 7 ') + 



1 - e 



[fi 2 + f 2 2 l 



± 4ea&xi£ 2 e 7 * cos(2cj£), 



(21) 



= ea&e" 7 ', 



(22) 



B(t) 



eb 2 



1 - e 



-27i 



(23) 



C(t) = 2e" 7t 



d> 2 



1 - e 



'1-e 



-7t\ 



1 - e 



(24) 



In the next section we will explore these solutions plotting the corresponding quantum trajectories for the generalized Werner 
state of Eq. ( Tl3T > in regions occurring entanglement or separability. 
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FIG. 1: Bohmian trajectories for the generalized Werner state of Eq. dl5t with a — b. For separable states with (a) e — 0.1 and (b) e = 1/3. 




FIG. 2: Bohmian trajectories for the generalized Werner state of Eq.d 151l with a = b. For entangled states with (a) e = 0.4 and (b) e = 1. 

IV. QUANTUM TRAJECTORIES FOR ENTANGLED STATES 

In this section we present our results regarding Bohmian trajectories for both separable and entangled states using the gener- 
alized Werner sta te g iven in Eq. ([151 , To quantify the entanglement present in this state we can either use the concurrence lfl6ll 
or the negativity lll7lfl8ll . thus we will use the concurrence as defined for two-qubit system: 

C( / 9)=max{0, v / A7- V^-V^- V 7 ^!}, (25) 

where Xk are the eigenvalues of the matrix pi2 (t) = cri<7^p* 2 (i)<T*cr5j arranged in decreasing order. 

As is well known, when disregarding losses and a = b, this state is separable for e = 1/3. In Figs. 1(a) and 1(b) we show the 
Bohmian trajectories for the generalized Werner state of Eq. ( fT3T > for a = b and e = 0.1, 1/3, corresponding to a separable state. 
In Figs. 2(a) and 2(b) we show the quantum trajectories when a = b and e = 0.4, 1, corresponding to an entangled state. From 
this sample of figures it can be seen that for separable state (e < 1/3) the amplitude of oscillations of each trajectory is smooth 
and relatively small as compared with the corresponding trajectory (same initial condition) for entangled states (e > 1/3), as 
advanced in Fig. 5(a) where we have plotted x(t) versus e for the same initial condition. From Figs. 2(a) and 2(b) we observe 
that, as the mixing parameter (e) increases, there is a corresponding squeezing of the trajectories in regions where both particles 
approach each other. This squeeze of trajectories, which clearly increases from Fig. 2(a) to Fig. 2(b), together with the increasing 
oscillation amplitude, works as a signature of entanglement for the generalized Werner state studied here. 

Using Eq. (l25l l. we can readily check that disregarding losses and considering a — 0.2, this state is separable for e < 0.56. 
In Figs. 3(a) and 3(b) we show the corresponding trajectories for separable states, while in Figs. 4(a) and 4(b) we show the 
trajectories for entangled states. Note that the same behavior as that for a = b = l/v2 can be observed in Fig. 5(b) that 
consider a = 0.2: given the same initial conditions, quantum trajectories for separable states are smooth and oscillate less than 
those corresponding to entangled states. For a complete mixture (e = 0), the trajectories are straight lines. 

Now let us see what happen to Bohmian trajectories in the presence of losses. Again using Eq. (l25l l. we can check that 
considering a = b and e = 0.4 this state is separable at jt — 0.15. Figs. 6(a,b) and 7(a,b) show the corresponding trajectories 



6 




FIG. 3: Bohmian trajectories for the generalized Werner state of Eq. H5\ with a = 0.2. Entanglement occurs for e > 0.56. (a) e = 0.1 and 
(b) e = 0.56. 




FIG. 4: Bohmian trajectories for the generalized Wemer state of Eq. d!5t with a — 0.2. Entanglement occurs for e > 0.56. (a) e = 0.7 (b) 
and e = 1.0. 

for damped states. As expected, the reservoir attenuates the oscillations, that turn to be more accentuated when the damping rate 
is larger. 

Note that these results for damped states show the same behavior as that for undamped states: given the same initial condi- 
tions, quantum trajectories for separable states are smooth and the amplitude of oscillations is less than those corresponding to 
entangled states, with the amplitude of oscillations going to zero when e — > 0. 

As a final remark, it is worth noting that Bohmian trajectories may provide a way towards detecting quantum separability 
of mixed quantum states. As a matter of fact, when the interpolator parameter e evolves from to 1, the trajectories changes 
from straight lines to curves with steep slopes. However, around the regions where the trajectories of the particles come closer 
together — thus interfering to a greater extent — their curvatures become very smooth. By focusing our attention on these regions 
of maximum interference between the trajectories, which takes place in Figs. 2 and 4 for Lot — (2n + l)%/2, n — 0, 1, 2, 
we thus observe that when the parameter e evolves from to 1, the slopes of the trajectories starts from 0, seems to reach 
a maximum value and then decreases due to the strong interference between the particles paths. We might suspect that the 
maximum curvature takes for the value of e that gives the separability condition for the density matrix. However, a problem arise 
when we set out to compute the curvature of the trajectories (in the specified regions) as a function of the parameter e, since the 
curvatures of the trajectories are different for different initial positions x a of the particles. 

V. CONCLUSIONS 

In this paper we have derived Bohmian trajectories for noninteracting bipartite states of damped harmonic oscillators under a 
thermal reservoir at finite temperature in a similar way to that of Vink's extension of Bell's beables lfl3ll . As an application, we 
have calculated the trajectories for a generalized Werner state dissipating at zero temperature in regions where the two systems are 
either entangled or separable according to Wooters' concurrence. Our results indicated that individual trajectories for entangled 
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FIG. 5: Bohmian trajectories for (a) a = b and (b) a = 0.2, for the same initial conditions and several values of e. The constant trajectory is 
for e = 0, while the oscilation amplitude is maximum for e = 1.0. The trajectories for (a) e = 1/3 and (b) e = 0.56 are indicated by a dashed 
lines. 




states differ slightly in the amplitude of oscillation as compared with those corresponding (same initial conditions) trajectories 
for disentangled states. We note, however, that according to our simulations, this difference is not enough to characterize 
unambiguously separability or entanglement, which is a global property of the system. This is so because the trajectories change 
continuously when the state changes from separable to nonseparable. We hope these preliminary results can encourage future 
research towards an eventual link between separability and Bohmian mechanics. 
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FIG. 7: Bohmian trajectories for the generalized Werner state of Eq. dl5b with a = 0.2, and w/7 = 10.0 for (a) for e = 0.7, where 
entanglement is suddenly lost at yt = 0.026 and (b) e = 1.0, where entanglement is lost at •yt = 0.23. 



